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At present there is no correct theory of evolution of ^2(0;, Q 2 ) at small x. 
It is a mixture of hard and soft pomeron exchange and perturbative QCD 
very successfully describes the evolution of the hard-pomeron component. 
This allows the gluon density to be calculated. It is somewhat different from 
what is conventionally supposed, but it leads to a clean PQCD description 
of the data for the charm structure function. Perturbative QCD breaks 
down for the evolution of the soft-pomeron component of F%(x, Q 2 )- 

1. Introduction 

The conventional treatment [2| of evolution expands the DGLAP split- 
ting matrix in powers of a s (Q 2 ). As we will explain, this is almost certainly 
wrong at small x, and at present we have no correct theory. However, when 
we combine PQCD with Regge theoryjH], this problem is partially solved and 
provides a very successful description of data, not only the complete proton 
structure function F2(x,Q 2 ) but also[I] its charm component F^x^Q 2 ). 

The proton's gluon density is larger at small x than is usually predicted, 
particularly at small Q 2 . A consequence of this is that PQCD evolution 
cleanly and successfully describes charm production at small Q 2 , even down 
to Q 2 = 0. See figure [TJ 

2. Regge theory the two pomerons 

At small x we make the fit 6 

F 2 (x, Q 2 ) = HQ 2 ) x~ e ° + h(Q 2 ) x~^ (1) 
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at each Q 2 for which there are data. We fix ei = 0.0808, the classical 
soft-pomeron value extracted from hadron-hadron scattering datajTJIH]- It 
turns out that, although the data are now highly accurate, they do not 
constrain the value of eo very closely. Good fits may be obtained with eo 
anywhere between 0.35 and 0.5. We call this the "hard-pomeron" term. 
While varying eo through its allowed range has little effect on the shape of 
the hard-pomeron coefficient function fo(Q 2 ), the large-Q 2 behaviour of the 
soft-pomeron coefficient function fi(Q 2 ) changes markedly; see figure^ 

For eo ~ 0.4 the data make /i(Q 2 ) go to a constant at large Q 2 . We 
assume[B] that fi(Q 2 ) has this behaviour and fit the available data for x < 
0.001. When we made the fit we used data from ZEUS 9 at small Q 2 
and from Hl^U] at larger Q 2 . There are now data at large Q 2 also from 
ZEUS [IT]. The best fit is now given by 

/o(Q2) = A) (i + QVQlY^i 2 MQ2) = Al WTWIW^ (2) 

with 

e = 0.4075 ei = 0.0808 

A) = 0.00227 Q = 2.88 GeV A 1 = 0.588 Qi = 768 MeV (3) 

See figure El We have already explained that the data do not constrain 
the value of eo very closely. We have given the parameters to this accuracy 
because their errors are strongly correlated. This set of values gives a x 2 
per data point significantly less than 1. 

It is an extremely economical fit: we included in it also data for pho- 
toproduction (figure below), which largely determine the value of A\, so 
that there are just 4 free parameters. If we multiply (Q) by (1 — x) 7 , which is 
a very crude way of ensuring that F2(x, Q 2 ) vanishes as x — > 1, and include 
a term corresponding to /2>«2 exchange, the fit agrees well with the data 
for larger x, even up to Q 2 = 5000 GeV 2 : figured 

If we try making a similar fit to the datajS] for the charm structure 
function, we find 6 that they correspond only to a hard-pomeron term. 
Further, the data are fitted well by assuming that the hard pomeron is 
flavour blind, so that for small x 

F^x,Q 2 ) = 0Af (Q 2 )x-^ (4) 

where f (Q 2 ) is defined in {TJ. The factor 0.4 is |/(| + | + | + |). Figure 
IS] shows 

V ' Q 2 2V ' x=Q*/(W 2 +Q 2 ) V ' 
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3. The DGLAP equation 

Define as usual the singlet parton densities 

*. ^ )=( :CE 4 ( (it, ^ ' , ) ' = MeVA 2 > (6) 
and take their Mellin transform with respect to x: 

u(N,Q 2 )= [ dxx N - l u(x,Q 2 ) (7) 
Jo 

Then the DGLAP equation reads 

^u(A,Q 2 ) = P(N,a s (Q 2 ))u(N,Q 2 ) (8) 

where P(iV, a s (Q 2 )) is the Mellin transform of the splitting matrix. 

The normal procedure is to expand P(iV, a s (Q 2 )) in powers of a s {Q 2 ). 
However, this is illegal when N is close to 0. This is well known. Compare, 
for example, the analogous expansion of the function 



= a s /2N -a 2 /8N 3 + ... (9) 

Although each term in the expansion is singular at A = 0, the function ip 
is not: the expansion is valid only for \N\ > a s . Similarly, the terms in the 
expansion of P(A, a s (Q 2 )) have singularities at N = which are surely not 
present in P(A, a s (Q 2 )) itself. Indeed, it is likely that P(A, a s (Q 2 )) has no 
relevant A-plane singularities at all. 

At any given Q 2 , expanding the splitting matrix in powers of the QCD 
coupling becomes invalid when one goes to sufficiently small x. 

At present, we have no other way to calculate. Luckily, if we introduce 
the two-pomeron parametrisation of the data we can partially rescue the 
situation. A fixed-power behaviour 

u(x,t)~ x - e (10) 

as occurs for each of the terms in (^Q) corresponds to 

If we insert this behaviour into the DGLAP equation (jHJ) and equate the 
coefficient of the pole at N = e on each side of the equation, we find an 
exact equation that describes how f(Q 2 ) evolves with Q 2 : 

^- t f(Q 2 )=P(N = e,a s (Q 2 ))f(Q 2 ) (12) 
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4. DGLAP evolution 

To calculate the evolution of the soft-pomeron term in we need 
P(N,a s (Q 2 )) at A = 0.0808. This is dangerously close to N = 0; we 
cannot make the expansion in powers of a s and do not know how to calcu- 
late the splitting matrix at this value of N . But for the hard-pomeron term 
we need P(iV, a s {Q 2 )) for N ~ 0.4, which is safely away from iV = and 
so the expansion should be valid. 

In order to solve the evolution equation (|12|) for the hard-pomeron coef- 
ficient functions f q (Q 2 ) and f g (Q 2 ), we chose Q 2 = 20 GeV 2 as our starting 
value. It does not matter what value we take, as long as it is not too small. 
We also assumed that, at this value of Q 2 , the conventional DGLAP anal- 
ysis of the data is correct for values of x down to about 0.01. That is, we 
assumed that the value of g(x = 0.01, Q 2 = 20) extracted from the data 
by MRSTfl] or CTEQ 2 is correct. Further because, as we have seen, the 
charm structure function is entirely hard-pomeron exchange at small x and 
because, as is well known, it is directly related to the gluon density, we de- 
duce that g(x, Q 2 ) at small x is entirely hard-pomeron exchange. Therefore 
we know the value of f g (Q 2 ) at Q 2 = 20. For the value of f q (Q 2 ) at Q 2 = 20 
we go to our fit to the data, that is we use (0) and Q. 

We then used 1)12(1 to evolve away from Q 2 = 20, in both directions. The 
result is rather astonishing: although the phenomenological function fo(Q 2 ) 
in (j2J) rises at large Q 2 as a power of Q 2 , while the solution to (|12|) rather 
rises as a power of logQ 2 , the two are in extraordinarily good numerical 
agreement over a wide range of Q 2 . This is shown in figure E3 We took 
4 flavours, with A LO = 140 Mev. We have found (3] that the output is the 
almost same whether we work to leading order in the coupling or next-to- 
leading. 

So, for Q 2 greater than about 5 GeV 2 perturbative QCD describes the 
evolution of the hard-pomeron component of F2(x, Q 2 ) extremely well. This 
is a significant success both for PQCD and for the two-pomeron description 
of F2(x,Q 2 ). It is no surprise that perturbative evolution breaks down at 
small Q 2 ; the DGLAP equation is supposed to be valid only for sufficiently 
large Q 2 . 



5. Gluon density 

According to what we have said, the proton's gluon density is 



xg(x,Q 2 ) = / s (Q 2 )x- S0 <Mx,Q 2 ) 



(13) 
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where 0(0, Q 2 ) = 1 and </>(x, Q 2 ) -> as x — > 1. A good numerical fit to the 
LO f g (Q 2 ) in the range 5 < Q 2 < 1000 is given by 

/ 8 (0 2 )^.32 (1+ ^ / ' 14)1+ , a/2 (14) 

At Q 2 = 20 GeV 2 the MRST or CTEQ LO gluon density is well described 
by 4>(x,Q 2 = 20) = (1 — x) 6 . Figure [7] compares our gluon distribution in 
LO with those of MRST and CTEQ at two values of Q 2 . The differences are 
evident and become even more pronounced in NLO, where our distribution 
is much the same but those of MRST and CTEQ are rather smaller and 
significantly less steep. 

In NLO, our gluon distribution turns out [3] to be almost the same as 
in LO. This is because we use the DGLAP splitting matrix only at N ~ 
0.4, where all but one of its elements are almost equal in LO and NLO. 
As x decreases, the conventional analysis involves the splitting matrix at 
progressively smaller values of N, where its elements are no longer the same 
in LO and NLO. At very small x the conventional analysis involves very 
small N, where it becomes illegal to use the perturbative expansion of the 
splitting matrix. 

At Q 2 = 200 GeV 2 and x = 0.0001 our gluon distribution is twice as 
large as that of MRST or CTEQ. The fact that our NLO gluon density is 
larger than the conventional one at small x will be significant for experiments 
at the LHC. 

The cleanest window on the gluon density will be provided by good 
data for the longitudinal structure function Fl(x,Q 2 ). Those data that 
exist depend on some assumed parametrisation to separate Fl from F2, for 
example reference |14j . Figure |H] shows that already at Q 2 = 20 GeV 2 there 
is a clear difference between our prediction and that of MRST. 

6. Charm production 

It is standard |15l 1161 117j that at small Q 2 the charm structure function 
F£(x, Q 2 ) should be calculated from photon-gluon fusion, 7*3 — > cc, to 
some fixed order in a s . This calculation introduces some assumed value for 
the charmed-quark mass m c . At large Q 2 a resummation to all orders in 
a s is needed, because of the presence of factors of powers of log Q 2 /m 2 . 
This resummation is achieved by changing at large Q 2 to the output from 
DGLAP evolution, where the charmed quark mass can now be neglected. 
The two calculations have to be matched at some value of Q 2 . The usual 
matching is done at a rather small value of Q 2 , of the order of m 2 , and is 
sensitive to exactly what value is chosen. 
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The thin lines in figureEJshow the result of the LO calculation of photon- 
gluon fusion with our gluon density and m c = 1.3 GeV. The results of an 
NLO calculation are almost the same, if we increase m c to 1.6 GeV. This 
calculation used the code of |18l I19j with our gluon distribution and m c as a 
free parameter. The thick lines in the figure are the fit @ which, as we have 
shown, agrees well with the output from DGLAP evolution for Q 2 greater 
than about 5 GeV 2 . So in our approach the two calculations match well 
over a range of Q 2 , from about 5 to 50 GeV 2 . 

7. The hard pomeron 

It used to be a central tenet of high energy physics that scattering am- 
plitudes are analytic functions of all their variables |2L)j. A consequence of 
this is that a singlarity that is present in an amplitude at large Q 2 survives 
when one goes to Q 2 = 0. In particular, if the hard pomeron is present at 
large Q 2 it should also be present at Q 2 = 0. 

This view is strongly reinforced by the charm-production data shown in 
figure [5J The W dependence at Q 2 = is the same as at higher Q 2 . The 
hard pomeron is not generated by PQCD evolution, though the evolution 
makes it relatively more important as Q 2 increases. Thus, in the fit we have 
described to F^ix, Q 2 ), there is a small hard-pomeron component already at 
small Q 2 , though its significance is masked by the much larger soft pomeron 
component. Both grow with increasing Q 2 , but the hard-pomeron compo- 
nent grows faster, and at each small value of x it dominates at sufficiently 
large Q 2 . 

Figure shows the data for the photoproduction cross section. The 
curve is our old fit[71|H], with no hard-pomeron term. The data are not 
yet good enough to test whether the fit is adequate or whether an extra 
component is needed such that indeed the hard pomeron is present already 
at Q 2 = 0. The same statement may be made of the LEP data |211 \2'2\ for 
cr 77 , which depend too heavily on Monte Carlo simulations that correct for 
poor acceptance to reach any conclusion. 

If the hard pomeron is present in the total cross section for photon col- 
lisons, is the same true for pp collisions? While probably the hard pomeron 
couples to a small object such as the photon with larger relative strength 
than to a large object such as the proton, there is some prospect that LHC 
data will show that there is a hard-pomeron component to a pp . 

What is the hard pomeron? Everybody agrees that the sharp rise in 
F2(x, Q 2 ) at small x discovered at HERA is a consequence of gluon exchange. 
Our own belief is that it is caused by glueball exchange and that the hard and 
soft pomerons are just glueball Regge trajectories. There is some evidence 
that this is true for the soft pomeron: there is a 2 ++ glueball candidate at 
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1926 MeV, exactly the right mass to be on the soft-pomeron trajectory 23 . 
Another 2 ++ glueball candidate [IS], a t 2350 MeV, could well be on the 
hard-pomeron trajectory |24|. 

At one time there was a hope that the power eo of 1/x, which is the hard- 
pomeron-exchange term, might be calculated from the BFKL and therefore 
that it is a perturbative effect. The soft-pomeron trajectory surely cannot be 
calculated from perturbative QCD. It may be that the glueballs on the hard- 
pomeron trajectory are heavy enough for their masses to be calculated from 
PQCD though, with the problems that have arisen with the BFKL equation, 
it is far from clear that PQCD can be used to calculate the intercept 1 + eo 
of the trajectory. 

8. Summary 

• The conventional approach to evolution needs modifying at small x 

• It can be corrected if we combine it with Regge theory 

• But only partly — we can only treat the hard-pomeron part 

• This is enough to extract the gluon distribution 

• The gluon distribution is larger at small x than has so far been supposed 

• It gives a good description of charm production 

• We want good data for the longitudinal structure function 



This research was supported in part by PPARC 
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Fig. 1. Data 5 for F|(x, Q 2 ) at Q 2 — 1.8 GeV 2 with theoretical curves from refer- 
ences £Q and 0] 
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Fig. 2. The hard and soft pomeron coefficient functions extracted from data 
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Fig. 3. Data from ZEUS^E] and HldHj with two-pomeron fit®. Q 2 ranges from 
0.045 to 45 GeV 2 . 




Fig. 4. The fit shown in figure 01 extended to larger x and Q 2 
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Fig. 5. Data[5] for a c {W) defined in (JSJ). The thick lines correspond to @, which 
coincides with the output from DGLAP evolution, and the thin lines are from 
photon-gluon fusion. 
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Fig. 6. LO and NLO calculations of the hard-pomeron coefficient function, together 
with the phenomenological fit J3J) 
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Fig. 7. The LO gluon density xg(x, Q 2 ) at Q 2 = 20 GeV 2 and 8 GeV 2 . The CTEQ 
and MRST curves are from the Durham Data Base|13j. 
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Fig. 8. The longitudinal structure function at Q 2 = 20 GeV 2 : data from Hl[T4] 
with our prediction (upper curve) and that of MRSTpQj 
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Fig. 9. 7p total cross section; the curve takes account of the exchange of the soft 
pomeron, / 2 and ai 



